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. , ABSTRACT 

' We present the results of a survey of radial velocities over a wide region extending from 

r ~ 10' out to r ~ 80' (~1.5 tidal radii) within the massive star cluster w Centauri. 
Q . The survey was performed with FLAMES@VLT, to study the velocity dispersion pro- 

|H ' file in the outer regions of this stellar system. We derived accurate radial velocities for 

C/5 I a sample of 2557 newly observed stars, identifying 318 bona-fide cluster red giants. 

■ Merging our data with those provided by Pancino et al. (2007), we assembled a final 

homogeneous sample of 946 cluster members that allowed us to trace the velocity dis- 
persion profile from the center out to r ~32 arcmin. The velocity dispersion appears 
to decrease monotonically over this range, from a central value of tr^ ~ 17.2 Km 
down to a minimum value of cr„ 5.2 Km s^^. The observed surface brightness pro- 

, file, rotation curve, velocity dispersion profile and ellipticity profile are simultaneously 

' well reproduced by a simple dynamical model in which mass follows light, within the 

I classical Newtonian theory of gravitation. The comparison with an N-body model of 

T^lj- i the evolution of a system mimicking ui Cen during the last 10 orbits into the Galactic 

' potential suggests that (a) the rotation of stars lying in the inner ~ 20' of the clusters 

, is not due to the effects of the tidal field of the Milky Way, as hypothized by other au- 

■ thors, and (b) the overall observational scenario is still compatible with the possibility 

that the outer regions of the cluster are subject to some tidal stirring. 
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1 INTRODUCTION 

u> Centauri is the most massive and luminous globular clus- 
ter (GO) of the Milky Way {M ~ 2.5 10** Mq; van de 
Ven et al. 2006). It is the only known Galactic GC which 
shows clear star-to-star variations in the abundance of iron- 
peak elements (Freeman & Rodgers 1975; Norris, Freeman 
& Mighell 1996) and indication exists for possible helium 
abundance variation among its stellar populations (Norris 
2004; Piotto et al. 2005). Together with a few other clusters 
(like M54 and NGC2419, for example), uj Cen fall well above 
the sharp upper envelope of the main distribution of GCs on 
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the size-luminosity plane (Mackey & van den Bergh 2005; 
Federici et al. 2007). These (and other) evidences suggest 
that u) Cen could be not a "genuine" GC but more likely 
the nuclear remnant of a dwarf galaxy that merged in the 
past with the Milky Way (Freeman 1993; see Bellazzini et 
al. 2008 for a discussion of the analogy with the case of M54 
within the disrupting Sgr dwarf galaxy). 

From the dynamical point of view, uj Cen is one of the 
most flattened among Galactic GCs (Meylan 1987) . The flat- 
tening is generally interpreted as due to the observed rota- 
tion, as the variations of ellipticity with distance from the 
cluster center correlates with the amplitude of the rotation 
curve (Merritt, Meylan, & Mayor 1997; Meylan & Mayor 
1986). Merritt et al.(1997) found that lo Cen can be well 
described as an isotropic oblate rotator in which the mass 
is distributed as the light. In a recent study Van de Ven 
et al. (2006) modelled (jj Cen with an axisymmetric imple- 
mentation of Schwarzschild's orbit superposition method. 
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They found that the system is close to isotropic inside a 
radius of about 10 arcmin and becomes increasingly tan- 
gentially anisotropic in the outer region, which displays sig- 
nificant mean rotation. These authors suggested that this 
phase-space structure could be caused by the effects of the 
tidal field of the Milky Way, as kinematical gradients are a 
typical outcome of tidal stresses (see Mufioz et al. 2008 and 
references therein). 

A number of studies on the kinematics of the cluster 
have been performed in the recent past (Meylan et al. 1995; 
Mayor et al. 1997; Norris et al. 1997; Reijns et al. 2006; Van 
de Ven et al. 2006), particularly focused on the innermost 
r < 20' region. Recently, Scarpa, Marconi & Gilmozzi (2003) 
measured the velocity dispersion profile of u Cen at larger 
distances (20' < r < 28') using UVES spectra of 75 candi- 
date cluster members. They found that the velocity disper- 
sion profile remains constant at large radii rather than de- 
crease monotonically, as expected for systems in which mass 
follows light. These authors suggested that this behaviour 
could be due to a breakdown of Newton's law in the weak ac- 
celeration regime, as proposed in the popular Modified New- 
tonian Dynamics framework (MOND, Milgrom 1983, 2008). 
However, it has been noted that actually this cannot be the 
case, as a flat velocity proflle at large radii is predicted by 
MOND only for clusters in the deepest MOND regime, i.e. 
where both the internal accelleration and that due to the 
gravitational field of the Milky Way are significantly lower 
that the MOND scale accelleration, ao ~ 1.2 x 10~^ cm/s^ 
(Baumgardt, Grebel & Kroupa 2005; Moffat & Toth 2008). 
On the contrary, at the position of u) Cen the external ac- 
celleration due to the Milky Way is > ao- In this regime, 
the predictions of MOND are similar to the Newtonian be- 
haviour (see also Moffat & Toth 2008; Haghi et al. 2009). 

Nevertheless, a flat velocity dispersion profile over a 
large radial range would be uncompatible with a mass- 
follow-light model (Gilmore et al. 2007) and may suggest the 
presence of a halo of Dark matter (DM) embedding the star 
cluster and driving the kinematics of the stars in the outer 
regions (Carraro & Lia 2000; Mashchenko & Sills 2005). In 
any case the available evidences are not sufficient to support 
the case of DM: as McLaughlin & Meylan (2003) demon- 
strated, Scarpa et al.'s data could be well described by sim- 
ple, self-consistent dynamical models without the need of 
MOND and/or dark matter. A larger kinematic survey in 
the region r > 20' is overdue to assess the behaviour of the 
dispersion profile in this range. 

A controversial issue is also represented by the possi- 
ble presence of signatures of tidal effects in the outskirts 
of this stellar system. If the present-day cluster is just the 
nuclear remnant of a larger system, it is reasonable to ex- 
pect some observable overdensity of stripped stars in the 
surroundings, if the latest phase of disruption occurred in 
the recent past (up to a few orbital period ago; see Combes, 
Leon & Meylan 1999; Dinescu (2002); Tsuchiya, Korchagin 
& Dinescu 2004; Mizutani, Chiba & Sakamoto 2003; Ideta 
& Makino 2004; Bekki & Freeman 2003). Leon, Meylan & 
Combes (2000) studied the 2-D structures of the distribution 
of stars around u Cen from 5.5° x 5.5° photographic films 
obtained with the ESO Schmidt telescope. This analysis ev- 
idenced the presence of a pair of tidal tails approximately 
oriented in the direction of the Galactic center, in projection. 
This result has been questioned by Law et al. (2003) who 



found that a) Leon's et al. tidal tails were strongly correlated 
with unhomogeneities in the reddening distribution and b) 
no extra-tidal component can be detected using the homoge- 
neous near-infrared photometry of the 2MASS survey. In a 
recent dedicated survey. Da Costa & Coleman (2008) derived 
the radial velocities (with an accuracy of ~ 10 Km s~^) of 
4105 stars located on a wide region out to a distance of 2 
tidal radii from the cluster center. They estimated that less 
than 0.7% of the total cluster mass is comprised between 
1 and 2 tidal radii, implying that the stripping process on 
the progenitor (if any) must have been largely completed at 
early epochs. 

In this paper we present the results of a survey of radial 
velocities in ui Cen performed with FLAMESQVLT, focused 
on the radial range 10' < r < 30' and aimed at the assess- 
ment of the kinematical properties in the outskirts of the 
cluster. The new dataset, combined with the homogeneopus 
sample of Pancino et al. (2007), allowed us to derive a ro- 
bust velocity dispersion profile from the center to a distance 
of 32' with an accuracy ^ 1 Km In §2 we describe the 
observations and the data reduction techniques. In §3 we 
present the metallicity and the velocity distributions. Sec- 
tion 4 is devoted to the description of the method used to 
derive the velocity dispersion profile of uj Cen. In §5 we com- 
pare the observed kinematics of this stellar system with a set 
of dynamical models and N-body simulations. Section 6 is 
devoted to the analysis of the sample of stars located in the 
outer region of uj Cen. Finally, we summarize and discuss 
our results in §7. 



2 OBSERVATIONS AND DATA REDUCTION 

The analysis presented here is based on two spectroscopic 
datasets: i) the mner sample: constituted of a sample of 
700 red giants selected from the photometry by Pancino et 
al. (2000) in the magnitude range 13 < B < 16 lying within 
15' from the cluster center and presented in Pancino et al. 
(2007; hereafter P07), and ti) the outer sample: constituted 
of a sample of 2557 red giants selected from the 2MASS 
catalog (Skrutskie et al. 2006) between 6.0 < K < 13.0, 
observed in 28 pointings at distances 10' < r < 80' from the 
cluster center (see Figure [l|. 

At r ~ 20' from the cluster center, the surface bright- 
ness is already fiy > 22.0 mag/arcsec^, implying that Red 
Giant Branch (RGB) stars of u Cen should be relatively 
rare and indeed not numerous enough to fill all the fibers of 
a single FLAMES plate. To achieve the maximum efficiency 
in selecting cluster members for the outer sample we gave 
the highest priority for fiber allocation to stars (a) lying 
within a selection box in color-magnitude diagram broadly 
enclosing the cluster RGB (see Fig. [5]), defined from the ac- 
curate infrared photometry by Sollima et al. (2004) , and (b) 
with UCAC2 (Zacharias et al. 2004) proper motions differ- 
ing by less than ± 8 mas yr~^ from the bulk proper motion 
of the cluster, as determined by van Leeuwen et al. (2000). 
These "a priori" criteria allowed to select some hundreds 
high-confidence target stars. The remaining fibers were po- 
sitioned on other stars selected from the 2MASS catalog in 
the same magnitude range, as shown in Fig. O below. 

Observations have been done with FLAMES (Pasquini 
et al. 2002) at the ESO VLT in Paranal, Chile, between 2003 
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May 22 and 28 {inner sample, P07) and between 2008 July 
13 and August 20 {outer sample). In both cases, FLAMES 
has been used in GIRAFFE mode, using the high-resolution 
(R ~ 22,500) grating HR13 (6120-6395 A) and reaching a 
signal-to- noise ratio of 50 — 300 pixel~^ , depending on the 
star magnitude. The inner sample has been already used by 
POT to study the rotation pattern of the stellar populations 
of ui Cen. The data reduction and calibration procedure is 
very similar to that adopted here for the outer sample (see 
below), and is described in detail in P07. 

The outer sample data have been reduced with the GI- 
RAFFE BLDRS (Base-Line Data Reduction Software) 4 
which includes cosmic-ray removal, bias subtraction, flat- 
field correction, wavelength calibration, and pixel resam- 
pling. Sixteen fibers have been dedicated to sky observa- 
tions in each exposure. An average sky spectrum has been 
obtained and subtracted from the object spectra by taking 
into account the different fiber transmission. The spectra 
have been then continuum-normalized and corrected for tel- 
luric absorption bands with IRAF. The spectra of three pro- 
gram stars are shown in Figure [2] to illustrate the quality of 
our data. 

Radial velocities have been derived through Fourier 
cross-correlation, using the fxcor task in the radial velocity 
IRAF package. The spectrum of each object has been cor- 
related with a high signal-to-noise template spectrum of the 
Geneva radial velocity standard star HD42807, retrived from 
the ELODIE archive (Moultaka et al. 2004). All spectra have 
been corrected for heliocentric velocitjQ. The average error 
on the derived radial velocity is cr„ ~ 0.5 Km s~^. Nine stars 
are in common between the inner and the outer samples: the 
average difference in radial velocity is A «r = —0.09 ± 0.16, 
fully consistent with no systematic shift between the two 
samples. Star by star comparisons with other datasets are 
also very satisfying. The average radial velocity difference 
of the 278 stars in common with Reijns et al. (2006) is 
A Vr = 1.48 ± 1.76, while for the 36 stars in common 
with Johnson et al. (2008) is A t;^ = 1.62 ± 1.47, i.e., 
consistent with null difference within the uncertainties (Re- 
jins et al. velocities were derived from spectra of resolution 
8500 < R < 17000; the resolution of Johnsons et al. spec- 
tra was R=13000; in both cases the typical uncertainty is 
- 1.5 Km s"^). 

In addition, metallicities have been also derived for the 
outer sample stars. A set of 15 iron lines has been selected 
from the database of Kurucz & Bell (1995) in the spectral 
range covered by our spectra. These lines have been iden- 
tified, whenever possible, on each spectrum and have been 
fitted with Gaussian functions. The integral of the differ- 
ence between the continuum and the line profile provided the 
equivalent width (EW) of each line. The abundance analysis 
has been performed using the latest version of the Kurucz 
(1979) model atmospheres and the MOOG line analysis code 
(Sneden 1973) to compute LTE abundances from individual 
EWs. The abundance has been derived from the model that 
best reproduced the observed EWs, for assumed values of 
temperature, gravity, and microturbulence velocity. Temper- 
atures have been derived from J-K colors adopting the color- 



The derived radial velocities of the outer sample stars are avail- 
able in electronic form at the CDS ( http:/ /cdsweb.u-strasbg.fr/| . 
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Figure 1. Map of the region sampled by the FLAMES obser- 
vations. North is up, east towards the right-hand side. The 28 
fields observed with FLAMES {outer sample) are shown as con- 
tinuous circles. The cluster center and tidal radius (from Trager 
et al. 1995) are indicated by the black cross and the dashed line, 
respectively. 



temperature transformations by Montegriffo et al. (1998), a 
reddening E{B -V) = 0.11 (Lub 2002) and the extinction 
coefficients by Savage & Mathis (1979). The gravity parame- 
ter log g has been then computed assuming a distance mod- 
ulus to UJ Cen of (m - M)o = 13.70 (Bellazzini et al. 2004; 
Del Principe et al. 2006). The microturbulence velocity has 
been initially set as 2 Km and then adjusted within a 
range of about 1 Km s~^ by minimizing the trend in the 
deduced abundances with EWs for Fe I lines. The average 
error on the derived metallicities is o"[Fe/£ri ~ 0.2 dex. The 
comparison with the metallicity determinations by Johnson 
et al. (2008) for the 36 stars in common indicates a difference 
of A[Fe/ H]our-j = —0.21 ± 0.25, indicating an acceptable 
agreement for the purposes of the present study, that is fo- 
cused on the kinematics in the outer regions of u} Cen. Here 
we derived metallicities of our stars to have an independent 
sanity check of the selection criteria adopted in §3 to identify 
cluster members. A full analysis of the abundance pattern 
of the selected stars will be the subject of a more detailed 
dedicated analysis, that will be presented in a forthcoming 
paper. 



3 RADIAL VELOCITY AND METALLICITY 
DISTRIBUTION 

In Figure |3] the radial velocities of the entire sample are 
plotted as a function of the distance from the cluster center. 
To calculate distances we deprojected right ascension and 
declination into X and Y coordinates with the following re- 
lations, suited for extended objects that are not close to the 
celestial equator (from van de Ven et al. 2006): 



4 A. Sollima et al. 




V, = 3.33 Km s ' 




[Fe/II] = -0.97 



Fe I Fe T 



Fe I Fe I 



K = B.166 vv = 230.57 Km s"' [Fe/H] = -0.B2 



K = 7.833 V, = 242.93 Km s-' [Fe/H] = -1.61 




Fe I Fe 1 



6320 6330 6340 6350 

A(A) 




r (arcmin) 



Figure 2. Spectra of the field star 2MASS13331099-4749262 (top 
panel), and of the uj Cen stars 2MASS13262366-4742424 (middle 
panel) and 2MASS13283381-4732055 (bottom panel). A sample 
of iron lines used in the analysis are also indicated. 



Figure 3. The radial velocities of the entire sample are plotted 
as a function of the distance from the cluster center. Grey points 
represent object rejected by the radial velocity selection criterion 
(i) defined in Sect. 14.11 Crosses represent object rejected by the 
other selection criteria. 



X — — ro cos 5osin{RA — RAo) 

Y = ro[sin S cos So — cos S sin So cos{RA — RAo)] 

where RAo and So are the coordinate of the cluster cen- 
ter and ro = IOSOO/tt is the scale factor to have X and Y in 
arcmin. At the adopted distance, 1 arcmin corresponds to 
1.6 pc. As can be noted, the clump of uj Cen stars is clearly 
visible at large positive velocities {vr > 190 Km s~^) while 
field stars are broadly distributed around the average veloc- 
ity < Vr >— —8.0 ± 30.2 Km s~^ , in good agreement with 
the predictions of the Galactic model by Robin et al. (2003; 
hereafter R03) for thin disc and thick disc stars in this di- 
rection, and very different from the typical velocity of u Cen 
(< Vr >= 233.2 ± 0.4 Km s"^). 

Figure H] shows the metallicity distribution obtained for 
the 338 stars belonging to the outer sample and with radial 
velocity 190 < Vr < 290 Km s~^ . A sharp and asymmet- 
ric peak at [Fe/H] ~ —1.8 can be seen, with a long tail 
extending toward higher metallicity {[Fe/H] > —1.4). Sec- 
ondary peaks are also visible at metallicity [Fe/H] ~ —1.3 
and [Fe/H] ~ —0.85, in agreement with the previous spec- 
troscopic determinations by Norris et al. (1996), Suntzeff & 
Kraft (1996) and Johnson et al. (2008). 
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Figure 4. Metallicity distribution of the 338 stars belonging to 
the outer sample with radial velocity 190 < Vr < 290 Km s~^. 
The distribution is obtained using a running box of 0.1 dex width 
and a step size of 0.05 dex. 



4 VELOCITY DISPERSION 
4.1 Selection criteria 



In order to construct the velocity dispersion profile of uj Cen 
we need to distinguish the bona-fide cluster stars from the 
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Galactic field contamination. In fact, although the field pop- 
ulation and LO Cen stars have very different bulk velocities, a 
small (but non neglibible) fraction of stars belonging to the 
Galactic halo could lie in the same velocity range of genuine 
cluster stars. This is particularly important in the external 
region where the cluster density drops quickly with radius, 
while the density of the field population remains constant 
(see p. 

To this aim, we excluded from the inner sample all the 
stars with radial velocity Vr < 190 Km (see P07) and 
adopted three different selection criteria for the outer sam- 
ple: 

• i) Radial velocity: The bulk radial velocity of lj Cen 
has been estimated to be < Vr >~ 233.2 ± 0.4 Km 
with a maximum velocity dispersion at the center of ~ 
20 Km s"^ (Meylan et al. 1995). Therefore, a first selection 
has been made by taking as possible members only stars 
with 190 < Vr < 290 Km Note that this selection 
is safe since the velocity dispersion is expected to decrease 
significantly at large distances from the center. Therefore 
this criterion is not expected to introduce any bias in the 
selected sample. 

• ii) Position in the CMD: Figure [S] shows the K vs (J — 
K) CMD of the outer sample stars. We considered as clus- 
ter members only the stars lying within the selection box 
enclosing the cluster RGB that we have originally adopted 
to chose high-priority targets, that is shown in Fig. [S] 

• Hi) Proper motions: We identified 245 stars in com- 
mon with the proper motion catalog by van Leeuwen et al. 
(2000). These authors provide for each star a membership 
probability calculated according to its position in the proper 
motions plane. On the basis of this additional information, 
we excluded all the 20 stars with a membership probability 
smaller than 50%. 

A total of 946 stars passed the above selection crite- 
ria, 628 belonging to the inner sample and 318 belonging 
to the outer sample. In Figure [3] the accepted and rejected 
stars are clearly indicated. The adopted approach is rather 
conservative as our aim is to obtain the most reliable sam- 
ple of cluster members to study the kinematics in the outer 
regions where the velocity dispersion is low and can be bi- 
ased by a limited number of spurious sources. On the other 
hand, while the second and, in particular, the third crite- 
rion listed above may exclude from the sample some genuine 
member, this cannot bias our results as these selections are 
independent of the radial velocity of the stars. Nevertheless, 
is interesting to note that using the R03 Galactic model as 
decribed in ^ below, we expect a very low degree of con- 
tamination. In particular, the predicted number of Galactic 
stars with 190 < Vr < 290 Km s^^ in our sample is 2.2±0.5 
for r 15', 1.4 ± 0.3 for 15' < r ^ 28', and 1.4 ± 0.3 in the 
range 28' < r ^ 44' (see [JB] for a discussion of the r > 44' 
region) . 

Another potential selection criterion is based on the 
metallicity. Note, however, that only Galactic halo stars 
are expected to contaminate the sample at velocities 



^ The choice of an asymmetric range has been performed to min- 
imize the field contamination in the region where the bulk of the 
field stars are located. 
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Figure 5. K vs. (J—K) CMD of tt) Cen. Grey points represent the 
stars rejected by the radial velocity selection criterion (j)defined 
in Sect. 14.11 Open points represent stars rejected by the other 
selection criteria. The adopted selection box is also shown. The 
clump of member stars around K ~ 11.7 is the RGB bump of the 
cluster (see SoUima et al. 2004). 

190 Km s~^ < Vr < 290 Km s~^. This Galactic compo- 
nent have a metallicity distribution ranging between —2.2 < 
[Fe/H] < -0.8 with a peak at [Fe/H] = -1.6 (Ivezic et 
al. 2008) Moreover, the contaminating halo stars lie in the 
foreground and background of uj Cen. Their different dis- 
tances translate into a wrong determination of their surface 
gravities and, consequently, spread out their original metal- 
licity distribution. As a consequence, the resulting metallic- 
ity distribution of the halo population largely overlaps the 
metallicity distribution of the u Cen stars. Therefore, this 
selection criterion is poorly efficient in discriminating the 
cluster membership. However, we considered also a sample 
retaining only the stars whose metallicity is close to the 
main peak of the metallicity distribution shown in Figure 
H (-1.95 < [Fe/H] < -1.55). In the following paragraph 
we calculated the velocity dispersion profile of lo Cen with 
and without adopting this last criterion, as a further check 
of our results. 

4.2 Method 

Once selected the bona-fide cluster member stars, as a first 
step we corrected for perspective rotation. In fact, since uj 
Cen has a large extension on the plane of the sky (with a 
diameter about twice that of the full moon) , its systemic mo- 
tion produces a non-negligible amount of apparent rotation 
(Feast 1961). We corrected for this effect using the following 
correction (see van de Ven et al. 2006): 

Avr = 0.001379 {X^x + Yhy)D Km s"^ 
to be subtracted to the observed Vr. i-ix = 3.88 and 
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fiY = —4.44 are the systemic proper motions in units of 
mas yr~^ (van Leeuwen 2000) and D = 5.5 is the distance 
to u) Cen in kpc (Bellazzini et aL 2004; Del Principe et aL 
2006). 

As a second step, we divided the sampled area in a set 
of concentric annuli at various distances from the cluster 
center, searching for the best compromise to have well pop- 
ulated and compact radial bins (except for the innermost 
r < 1' bin that is not relevant in this context). Within each 
annulus we calculated the value of that maximize the 
log-likelihood 



Table 1. Velocity dispersion profile of w Cen. 
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where Vi and ai are the radial velocity and its corresponding 
error for the N stars contained in the annulus and < v > 
the average cluster velocity (Walker et al. 2006). As a fur- 
ther selection, in each annulus we excluded the outliers stars 
whose radial velocity differs more than 3 times the local 
from the cluster bulk velocity (see Figure [6l Bellazzini et al. 
2008). 

The derived velocity dispersion profile is listed in Table 
1 and shown in Figure [7] for the case of the sample selected 
with (filled points) and without (open points) the metallic- 
ity selection criterion. As can be noted, no significant differ- 
ences are evident between the two profiles. For this reason, 
in the following we refer always to the sample not selected in 
metallicity. The obtained profile shows a central velocity dis- 
persion of av = n in agreement with the pre- 
vious estimates by Meylan et al. (1995), Mayor et al. (1997) 
Norris et al. (1997) and van de Ven et al.(2006). The mea- 
sured velocity dispersion appears to decrease monotonically 
with distance out to r ~ 26' from the cluster center, where 
it reaches a minimum of (j„ — 5.2tg gg7fm . A small, non 
statistically significant, increase of the velocity dispersion is 
noticeable in the last bin, at a distance of ~ 32'. 

In Figure |5] the obtained velocity dispersion profile is 
compared with the profiles obtained by van de Ven et al. 
(2006) and Scarpa et al. (2003). Note that our profile agree 
very well with that obtained by van de Ven et al. (2006) in 
the overlapping region r ^ 20'. At distances r > 24' the be- 
haviour of our profile differs from that obtained by Scarpa 
et al. (2003). In particular, while the profile by Scarpa et al. 
(2003) fiattens around a value of ct^ ~ 8 Km s~^, our profile 
continues to decrease. The statistical disagreement between 
the two profiles is ~ 2a, hence it cannot be considered as 
statistically significant. The number of target stars in range 
20' r < 28' and the accuracy of the radial velocity esti- 
mates of the present study is similar to Scarpa et al. (2003) . 
However, Scarpa et al. (2003) do not report (a) if their veloc- 
ities have been corrected for perspective rotation (§4.2) and 
(b) the details of the adopted selection criteria; in particular 
it is not said if the selection in Vr takes into account the local 
value of the velocity dispersion at the considered radius. In 
any case, our results indicate that the fiattening of the veloc- 
ity dispersion curve shown by Scarpa et al. (2003) may not 
refiect a real kinematic feature of the cluster but a chance 
fiuctuation, instead. The (doubtful) increase of the disper- 
sion in the < r >~ 32' bin observed here, if real, would be 
more compatible with the onset of tidal heating that with 
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Figure 6. Radial velocity distribution as a function of the 
distance from the cluster center for stars with 190 < Vr < 
290 Km . The adopted bins of variable size are indicated to- 
gether with the local 3cr range. Grey points represent the object 
rejected by the selection criteria defined in Sect. 14.11 



the effects of MOND or DM (see see also Munoz et 

al. 2008, and references therein). Note that this last bin is 
located outside the radial range covered by Scarpa et al.'s 
data. 



5 COMPARISON WITH THEORETICAL 
MODELS 

To check if canonical dynamical models are able to repro- 
duce the observed kinematics of uj Cen, in the following sec- 
tions we compare our data with a set of analytical mod- 
els (Wilson 1975). N-body simulations have also been per- 
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Figure 7. Velocity dispersion profile of u) Cen. Open and filled 
points refers to the sample selected and not selected in metallicity, 
respectively. 



20 

r (arcmin) 



Figure 8. Comparison between the velocity dispersion profile 
presented in this paper (filled points) with those by van de Ven et 
al. (2006; open points) and by Scarpa et al. (2003; open triangles). 



5.1 Analytical models 

We attempt to reproduce our observations using Wilson 
(1975) models. These models have the advantage of taking 
into account the cluster rotation, predicting simultaneously 
the surface brightness, the velocity dispersion, the rotation, 
and the ellipticity profiles. These models adopt the distribu- 
tion function 



f{E, J) = 







1+E)e 



a E <o 



formed to study the effects of the tidal field of the Milky 
Way on different models of the cluster. 



Any choice of the a-dimensional central potential Uo, 
and of the two free parameters of the above equation (/3 and 
(), yelds a different rotating model in virial equilibrium. The 
adoption of a core radius Vc and of the inclination angle re- 
spect to the plane of the sky i provide an actual realization 
of the model that can be compared with observations. We 
explored the parameter space to find out one set reproduc- 
ing simultaneously i) the surface brightness profile, ii) the 
velocity dispersion profile, in) the rotation curve and iv) 
the ellipticity profile, giving lower weight to this last item. 
We adopt the velocity dispersion profile derived in Sect. 14.21 
the surface brightness profile by Trager, King & Djorgovski 
(1995) and the ellipticity profile by Geyer, Nelles & Hopp 
(1983). Since the angle of rotation in the plane of the sky 
of ui Cen is ~ 0°, the rotation curve corresponds to the 
X — Vr diagram. 

The best-fit model has been found by fixing the incli- 
nation angle i = 48° (van de Ven et al. 2006) and assuming 
Uo = -6.3, f3 = 0.8, C = 0.4 and = 2.8'. The ratio be- 
tween the rotational kinetic energy and the absolute value 
of the potential energy of this model is T/|W^| — 0.09, indi- 
cating that it is stable against non-axisymmetric perturba- 
tions (Ostriker & Peebles 1973). The fit of the observational 
data with this model is shown in Figure |9] As can be noted, 
the agreement is good in all the four diagrams. In partic- 
ular, the model well reproduces all the observed kinematic 
properties of the cluster (rotation and dispersion) over the 
entire radial range covered by our datctfj. This means that 
the velocity dispersion and rotation profiles of u) Cen are 
fully consistent with an equilibrium model in which the mass 
distribution follows the light, in the framework of classical 
Newtonian gravitation. This is in good agreement with the 
results by McLaughlin & Meylan (2003), that were limited 
to the velocity dispersion profile. 



5.2 N-body simulations 

As described in Sect. [TJ there is no strong evidence for the 
presence of a significant population of tidally stripped stars 
around uj Cen. The most thorough study in this sense is the 
one by Da Costa & Coleman (2008) that, from their exten- 
sive radial velocity survey, concluded that less than 0.7% 



^ The marginal disagreement between the predicted and observed 
velocity dispersion profile in the innermost region could be in- 
duced by many second-order effects in modelling the shape of 
the velocity distribution in the internal cluster region (e.g. non- 
rotational anisotropy) and/or by the possible presence of an 
intermediate-mass black hole object in the center of u) Cen (as 
proposed by Noyola, Gebhardt & Bergmann 2008) 
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Figure 9. Comparison between the best-fit Wilson (1975) model (solid lines) and a) the projected surface brightness (from Trager et 
al. 1995) b) the cllipticity profile (from Geyer et al. 1983) c) the rotation curve and d) the velocity dispersion profile presented in this 
paper. In panel c) the la,2a and 3(T contours predicted by the model are also shown as dashed lines. 



of the total cluster mass is comprised between 1 and 2 tidal 
radii. This result imply that the stripping process must have 
been largely completed at early epochs. Therefore, if uo Con 
is indeed the remnant of a larger system destroyed long time 
ago, it must have lost a large part of its stars in a orbit with 
a larger apogalacticon with respect to the present day, be- 
cause of orbital decay due to dynamical friction (Zhao 2002; 
Bekki & Freeman 2003; sec also Tsuchiya et al. 2004; Mizu- 
tani et al. 2003 and Ideta & Makino 2004). This implies that 
the present-day system should be approximately in equilib- 
rium with the tidal field of the Milky Way and lie on a stable 
orbit, as its mass is now sufficently low to make the effect of 
dinamical friction negligible. 

We adopt this view as a working hypothesis for a simple 
N-body simulation aimed at studying the spatial distribu- 
tion and the kinematical properties of the stars eventually 
stripped from the cluster by the tidal field of the Galaxy 
while the system is in the present-day stable orbit (see Ideta 
& Makino 2004). The very good knowledge of the distance 
and 3-D motion of uj Cen (van Leeuwen 2000) gives a high 
predictive power, under the limitations of the assumed hy- 
pothesis. For more detailed and extensive studies see Bekki 
& Freeman (2003), Tsuchiya et al. (2004), Ideta & Makino 
(2004), and references therein. 

The simulation was performed with falcON, a fast 



and momentum-conserving tree-code (Dehnen 2000, 2002), 
within the NEMO environment (Tauben 1995). Gravity 
was softened with the kernel 'P2' of Dehnen (2000), with 
a softening length of 2 pc. The model cluster was let to 
evolve within the three-component (bulge-f disc-f halo) static 
Galactic potential 2b of Dehnen & Binney (1998), which has 
VrotiRo) ~ ^^-'^ * total mass within 100 kpc 

of ~ 1 X 10^^ M0. Lj Cen was represented by 50000 particles 
distributed according to an equilibrium King (1966) model 
having Wq = 6.0, n = 100 pc and M = 5 x 10^ Mq. This 
model has concentration (c) and tidal radius (rt) values sim- 
ilar to the real cluster; the central velocity dispersion is also 
very similar, a ~ 17.0 Km , but it should be recalled 
that the initial velocity distribution is completely isotropic, 
lacking any intrinsic rotation, f alcON is best suited to follow 
the evolution of non-coUisional systems. While the effects of 
two-body encounters have been shown to be negligible in 
u> Cen (Ferraro et al. 2006; Sollima et al. 2007), this may 
not be the case in the innermost region of our model. To 
minimize undesired effects due to two-body encounters in 
the densest part of the model we have adopted a soften- 
ing length approximately as large as the cluster core radius. 
For this reason the evolution in the innermost few Tc of the 
model should not be considered as a good representation 
of the evolution of the real cluster. On the other hand our 
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model is fully adequate to study the reaction of the outer 
regions to the tidal strain of the Galactic potential, i.e. the 
morphology and the kinematics of possible tidal tails, that 
is what we are interested into. Note also that the 2-body re- 
laxation time at the half-mass radius of our model (Binney 
& Tremaine 1998, their Eq. 4-9) is a factor of ~ 4 larger 
than the duration of our simulation, hence, in any case, the 
effects of particle encounters on the global properties of the 
model should be negligible. 

First of all we integrated backward in time, within the 
adopted potential, the orbit of a test particle whose initial 
position and 3-D velocity coincide with those of uj Cen at 
the present epoch (position from Bellazzini et al. 2004 and 
velocity from van Leeuwen et al. 2000). The orbit has a pe- 
riod of ~ 80 Myr, pericentric and apocentric distances of 
1.5 kpc and 6.5 kpc, respectively, very similar to the or- 
bits derived by Mizutani et al. (2003) and Ideta & Makino 
(2004). We launched our N=50000 particle model from the 
apocenter at the epoch t=0.972 Gyr in the past. With this 
Eissumption the cluster evolves within the adopted potential 
for more than 10 orbits from that time to the present epoch 
(t=0). At the end of the simulation the model is within 40 
pc of the current position of u) Cen, and the 3-D velocity 
is within 5% of the real cluster velocity. The distribution 
of the particles at the end of the simulation (see Fig. llOf) 
is in good agreement with the results by Ideta & Makino 
(2004) and Mizutani et al. (2003). The remnant have devel- 
oped two tidal tails extending over the whole orbital path. 
However, the tidal radius of the bound remnant is only 10% 
smaller than the initial model and there is a mere 8% of 
the particles that lie outside that radius after 10 perigalac- 
tic passages. This confirms that uj Cen is currently nearly at 
equilibrium with the tidal field of the Milky Way and that 
any tidal tail originating from the present day cluster should 
be very weak and hard to detect. In particular, it is interest- 
ing to note that only 0.4% of the particles are found in the 
annulus between 1 rt and 2 rt, in excellent agreement with 
the results of Da Costa & Coleman (2008). Rescaling to the 
observed total luminosity of the cluster the average surface 
brightness in that annulus due to tidally stripped stars is 
Hv ~ 30.0 mag/arcsec^. 

In Figure [TT] the final surface brightness and the veloc- 
ity dispersion profiles of the simulated satellite are shown. 
As can be noted, the surface brightness profile deviates from 
the King profile at r > 40', following the power-law trend 
typical of tidally stripped stars (see Johnston, Sigurdsson & 
Hernquist 1999; Mufioz et al. 2008, and references therein). 
This deviation occurs at a very low surface density level 
(about 3.5 orders of magnitude lower than the center), just 
outside the limit reached so far by any photometric inves- 
tigation. Note that neither King models nor Wilson's ones 
can reproduce this external behaviour. 

Also the velocity dispersion profile shows a departure 
from the predicted King profile. However, in this case the 
drift already occurs at r ~ 25' yielding a fiatter profile at 
larger distances. In any case, the derived profile is compati- 
ble with the observed one , within the uncertainties. Hence, 
the available observations appear still compatible with the 
tidal stripping scenario described by our N-body model. 
Fig. [TT] suggests that the "easiest" way to look for the ob- 
servational signatures of the ongoing formation of tidal tails 
would be to reliably extend the surface brightness profile 
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Figure 10. Position of the particles at the end of the N-body 
simulation A (present epoch) in galactocentric coordinates. The 
X-Y plane coincides with the Galactic Plane, the X-Z plane is 
perpendicular to it and passes through the Sun. Note the excellent 
agreement with the much finer simulation by Ideta &; Makino 
(2004; the upper right panel of their Fig. 2). 



using abundant and easy-to-select tracers as Main Sequence 
stars, searching for the power-law branch of the profile. 



6 THE OUTER REGION 

In Sect. |4] we presented the velocity dispersion profile of cu 
Cen calculated out to a distance of ~32 arcmin, and com- 
prising stars out to r = 44' from the cluster center. A more 
careful analysis is needed when considering the region at 
r > 44' i.e. close to the cluster tidal radius. In this re- 
gion, only seven stars passed the selection criteria defined 
in Sect. 14. Il The resulting velocity dispersion of this sample 
is ay = 36.0l7^5^ -^"^ about a factor of 2 larger than 
the central value, and a factor of 5 larger than the outermost 
point of the profile, at < r >~ 32. It is worth noticing that 
although the adopted selection criteria are efficient in elim- 
inating the contaminating field stars in the inner region of 
the system, their efficiency decrease at larger distances, as 
the surface density of uj Cen stars plunges below the back- 
ground level due to Galactic field stars. 

To check the membership of these sample of stars we 
performed various tests. First, we compared our sample with 
the predictions of the R03 Galactic model. Although this 
model is only a statistical representation of the Galactic 
field population (see Drimmmel et al. 2005), it gives a good 
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Figure 11. Surface brightness (left panel) and velocity dispersion (right panel) profiles for the bound remnant at the end of the N-body 
simulation A (filled circles). The dashed lines are from the King (1966) model that best fits the final surface brightness profile. 



estimate of the order of magnitude of the field contamina- 
tion. To minimize the effect due to smaU number statistics, 
we retrived (and merged together) two simulated catalogue 
covering an area of 5 sqare degrees each, in the direction 
of OJ Cen. The Vr distributions of simulated stars approxi- 
mately lying in the same magnitude and colour range of the 
observed stars is well fitted by a sum of two Gaussian dis- 
tributions, one for the main peak of thin disc -I- thick disc 
Galactic stars at —150 < Vr < 150 Km s~^, the other for 
the Halo population. Once fixed the density normalization 
of the two components, the number of stars expected in our 
sample, in any given range of r and is easily obtained by 
re-normalizing the Gaussian of the main component to the 
number of observed stars having —150 < Vr < 150 Km 
in that range. Adopting this procedure, the expected num- 
ber of field stars with 190 < Vr < 290 Km at distances 
r > 44' turns out to be Nfieid = 4.0 ± 0.9. Therefore, ac- 
cording to this test, about half of the outermost stars are 
likely non-members. Galactic Halo stars. 

Another possibility is to estimate the number of cluster 
stars expected at such distance by means of the cluster den- 
sity profile. Indeed, in a given radial bin (defined between 
ri and r2) the expected ratio between the number of cluster 
stars (Nmem) and field stars (Nfi^id) can be written as 

M 2TTrpmemdr rpmemdr 
^^mem J r\ ' ^ Jri ' 

Nfieid J^^ 2nrpfieiddr J^^ r dr 

where Pmem and p field represent the projected surface 
density of member stars and field stars, respectively. In the 
above equation the density of field stars has been assumed to 
be constant over the field of view of the cluster. We counted 
Nfieid as the number of stars with 190 Km < Vr < 
290 Km in each radial bin and estimated the constant 
C by best-fitting the ratio Nmem /Nfieid measured in the in- 
nermost bins (at r < 44') with the prediction of the best-fit 
Wilson (1975) model described in Sect. 15.11 Then, we es- 
timated the number of cluster members stars in the most 
external bin by using eq. [T] and Nfieid measured in this 
bin. Following this procedure we estimated Nmem = 0.4 at 



r > 44'. However, as shown in Sect. 15.21 the interaction be- 
tween u) Cen and the Milky Way would produce a departure 
from the King profile in the external region of the cluster, 
implying a higher surface density. Therefore, the number of 
cluster members stars in the external region is expected to 
be larger than this estimate. We performed the above pro- 
cedure also assuming at r > 40' a power-law density profile 
with index a — —3.8 (i.e. the slope predicted by our N-body 
simulation). Also in this case, the expected number of clus- 
ter members results Nmem = 1.5. Thus, according to this 
test, the majority of the high- velocity external stars should 
be field stars. 

It is interesting to note that none of the seven consid- 
ered stars pass the metallicity selection criterion defined in 
Sect. 14.11 A Kolmogorov-Smirnov test gives a probability 
smaller than 0.2% that the metallicity of these stars can be 
randomly extracted from the metallicity distribution of the 
entire sample of cluster members stars shown in Figure |4] . 
Given all the result above we can safely conclude that most 
of (if not all) the seven stars passing our selection criteria 
at r > 44' are likely not member of the cluster. 

Having established that, we draw the attention of the 
reader on a curious occurrence. In Figure [T^ the radial ve- 
locities of the stars with 190 < Vr < 290 Km is shown 
as a function of the right ascension for the observed sam- 
ple and for the particles of the N-body model described in 
Sect. 15.21 Note that the most external stars are not evenly 
distributed in the considered plane. Most of the r > 44' 
stars lying to the West of the cluster have Vr larger that 
the systemic cluster velocity, while the opposite is true for 
r > 44' Eastern stars. This apparent trend of mean Vr with 
RA is in the same sense as that expected to be imprinted 
by the Galactic tidal field on stripped stars, clearly visible 
in the N-body model (lower panel of Fig. I12|l . By means of 
Monte-Carlo extractions from the synthetic R03 catalogue 
described above, we estimated that the chance occurrence 
of such an asimmetry in Halo stars is lower than 5%. If the 
observed gradient would be revealed as real using a much 
larger sample of stars in this range of distances, this would 
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Figure 12. The radial velocities of the stars with 190 < Dr < 
290 Km is shown as a function of the right ascension for 
the observed sample (top panel) and for the N-body simulation 
described in Sect. |5^ (bottom panel). 

make very difficult to interpret most of these stars as be- 
longing to the Galactic Halo and the hypothesis of a tidal 
tail must be reconsidered. Another interesting issue emerges 
from this comparison: in the innermost region of the cluster 
(201.0° < RA< 202.3° in Fi ifT2)) . the clear S-shaped trend 
due to the cluster rotation is clearly visible In particular, in 
this case the Eastern branch of the rotation curve has larger 
average Vr with respect to the Western branch. This trend 
is just in the opposite direction of what produced in the N- 
body simulation by the Galactic tidal strain. This strongly 
suggest that the rotation in the inner 30' of lo Cen is not 
induced the interaction with the Milky Way, as previously 
suggested by van de Ven et al. (2006), but it is an intrinsic 
property of the system, instead. 



consistent with those obtained by previous authors within 
r ~ 20' (Meylan et al. 1995; Mayor et al. 1997; van de Ven 
et al. 2006) but differs from that derived by Scarpa et al. 
(2003) in the most external region of the cluster. In partic- 
ular, at odds with the Scarpa et al.'s analysis, we did not 
find any sign of flattening of the velocity dispersion profile 
in the region between 20' < r < 28'. 

The observed velocity dispersion profile of u Cen, as 
well as its structural and kinematical properties is well repro- 
duced by a Wilson (1975) equilibrium model in which mass 
follows light. Therefore, with the present data we find no 
evidences of neither the presence of dark matter nor MOND 
effects, at least in the inner 30' from the center. 

The majority of the seven candidate members we found 
at r > 44' are consistent with being iterlopers from the 
Galactic Halo. The only piece of evidence that is in marginal 
contrast with this view is that they seem to display a veloc- 
ity gradient similar to that expected for star tidally stripped 
from LU Cen, according to the prediction of a simple N-body 
model. The comparison with the same N-body model sug- 
gest that the well-known rotation pattern observed in the 
inner 30' of u Cen is not produced by the Galactic tidal 
strain, but it should be intrinsical to the stellar system, in- 
stead. 

The present study, mainly focused on the kinematics in 
the range 20' < r < 30', as well as that performed by Da 
Costa and Coleman (2008), focused on the search of extra- 
tidal stars (60' < r < 120'), reached the limits that can be 
achieved using RGB stars as tracers. Future studies aimed 
at assessing the shape of the velocity dispersion profile in the 
range 30' < r < 60' and trying to identify extra-tidal stars 
must rely on much more abundant Main Sequence stars. 
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7 SUMMARY AND CONCLUSIONS 

In this paper we present high accuracy radial velocities for 
a sample of more than 2500 stars in the direction of u> Cen. 
Adding the homogeneous sample by P07 and applying vari- 
ous un-biased selection criteria, we selected a sample of 946 
bona-fide cluster's members displaced over a wide region up 
to a distance of ~ 80' from the cluster center. 

We derived a very reliable velocity dispersion profile in 
the range 1.5' ^ r 28'. In this range the velocity dis- 
persion appears to decrease monotonically from a central 
value of a„ ~ 17.2 Km down to a minimum value of 
~ 5.2 Km s~^. A non-statistically significant raise to 
~ 7.0 ± 1.0 Km is observed in an outer extended ra- 
dial bin (28' ^ r < 44'), that must be considered with great 
caution because of the uneven radial distribution of the 26 
candidate member stars included. The obtained profile is 
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